Neutrinos from Pion Decay at Rest to Probe the Proton Strangeness 
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The study of the neutral current elastic scattering of neutrinos on protons at lower energies can be 
used as a compelling probe to improve our knowledge of the strangeness of the proton. We consider 
a neutrino beam generated from pion decay at rest, as provided by a cyclotron or a spallation 
neutron source and a 1 kton scintillating detector with a potential similar to the Borexino detector. 
Despite several backgrounds from solar and radioactive sources it is possible to estimate two optimal 
energy windows for the analysis, one between 0.65 — 1.1 MeV and another between 1.73 — 2.5 MeV. 
The expected number of neutral current events in these two regions, for an exposure of 1 year, are 
enough to obtain an error on the strange axial-charge of 6 x 10~^, more than 10 times smaller than 
available at present. 

PACS numbers: 25.30.Pt; neutrino scattering. 14.20.Dh; properties of protons. 



Introduction: The neutrino-proton elastic scatter- 
ing was proposed by Weinberg fT] as a tool to inves- 
tigate the neutral currents (NC). However, this reac- 
tion depends upon 'proton strangeness'~i.e., the strange- 
quark contribution to the axial form-factor of the proton- 
a quantity that at present is not reliably known. Sev- 
eral theoretical as well experimental investigations of 
neutrino-proton elastic scattering [SHl] tried to probe this 
quantity. The best experimental result [3] found a value 
for the proton strangeness similar to the one predicted in 
[2], but the 90% C.L. range is compatible with zero. 

Proton strangeness matters in many contexts, e.g.: for 
the measurement of supernova neutrinos f5] , for the cou- 
plings of certain dark matter candidates to nucleons [7] , 
for the polarized parton distribution functions, e.g., [S]. 
However, the flavor singlet term measured by deep inelas- 
tic scattering of charged leptons includes an anomalous 
gluonic contribution [5]; unless one postulates that this 
contribution is zero [lOj . a considerable amount of addi- 
tional labor is needed to extract proton strangeness, see 
6-g-, By contrast, a measurement of elastic neutrino- 
proton scattering probes this quantity directly. 

In the recent years, neutrino detectors underwent ma- 
jor developments. The Borexino detector, in particular, 
has reached an unprecedented low energy threshold for 
the detection of solar neutrinos ^12j, implying the possi- 
bility to use highly pure scintillator detectors to measure 
protons with kinetic energies as low as 1.3 MeV. Thus, 
these technological achievements pave the way for new 
approaches to measure precisely the elastic scattering of 
neutrinos onto protons. 

We study the possibility to determine the proton 
strangeness measuring low energy neutrino-proton scat- 
tering. The advantages of the low energy regime are dis- 
cussed in the next section. An artificial neutrino beam 
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produced by pion decay at rest, as in the cyclotrons con- 
sidered by the DAE5ALUS proposal [T3J [H] or in the 
SNS facility [T5] is considered. Recall that the neutrino 
spectra from pion decay are very well known, and have 
already proved being useful for measurements of other 
cross sections [TB]. The events are observed by a scin- 
tillating detector. We consider a 1 kton mass detector, 
located at 100 m from the source and with the perfor- 
mances of Borexino 12 . We explore the potential of this 
experimental setup and determine which sensitivity can 
be reached with 1 year of exposure. We discuss the de- 
pendence of the sensitivity on the experimental features, 
showing how to rescale our findings for different config- 
urations; i.e., by varying the detector mass, the distance 
to the source, the time of data taking, etc.. 

Advantages of Low Energies: The neutral current 
elastic scattering transition i/(fc) + p{p) v{k') + p{p') 
is described by the invariant matrix element 

G 

^ u'^i'^PlUu u'pTa{q)up (1) 

with 

Ta{q) = laFiiQ^) -f i^F^iQ^) - 7a75G.(Q2) (2) 
imp 

and similarly for antineutrinos, where q ^ p — p' la the 
transfered momentum and = —q^. The vector form 
factors Fi,2((5^) can be obtained by applying the con- 
served vector current hypothesis neglecting the strange 
vector form factors. The most important quantity for us 
is the axial form factor Ga(Q^). This can be determined 
from the analogous term in the charged current interac- 
tions, but only up to the SU(2) isosinglet term, due to the 
strange quark axial current. Its effect has been described 
by the parameter 77: 

^ ,^2^ _ 9M X (1 + 7?) , . 

where the dipole parameterization of the axial form factor 
is used and we have (7a(0) = 1.267. Adopting the usual 
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FIG. 1. Variation of TZ{E^) with the neutrino /antineutrino 
energy E^. The black solid line corresponds to r] = 0.12, while 
in the red dashed and brown dotted line r) has been set to zero. 
The red dashed line represents the effect of a relatively higher 
axial mass Ma = 1.2 GeV in the numerator of Eq^ See the 
text for detail. 



definition {p\qj°'j^q\p) = Upj^-j^Up ■ Aq the electroweak 
standard model predicts Ga(0) = Au/2 - Ad/2 - As/2, 
which agrees with Eq. Q setting 5a (0) = Am — Ad and 

r] = -As/g^iO). 



In this work, the term 'proton strangeness' indicates rj; 
the conversion to As is evidently trivial. As recalled in 
the introduction, the BNL 734 experiment [3] obtained 



T] = 0.12 ±0.07 



(5) 



that implies 77 = — 0.25 at 90% C.L., compatible with 
zero. Moreover, due to the use of neutrinos with ener- 
gies around 1 GeV, the experimental determination of 77 
is strongly entangled with the high behavior of the 
axial form factor, that is to some extent uncertain. Con- 
versely, the impact of the 'axial mass' on the cross section 
is greatly reduced if lower energies are considered. 

This is clear from Fig.jl] where we display the following 
cross section ratio: 



n{E,) 



y } ~r (7 up 



(6) 



In the denominator of Eq. we put the reference cross 
section, as given in with 77 = and Ma = 1.032 
GeV. In the numerator we modified it as follows: 1) In 
the black curve, we use the central value rj = 0.12 of 
Eq. ([5]) instead. Evidently, this causes an increase in the 
cross section. 2) In the red (dashed) curve, we set 77 = 
and use Ma = 1.2 GeV. By increasing Ma, the axial 
form factor and thus the cross section increase. 3) In the 
brown (dotted) curve, we set 77 = and omit the vector 
current contribution altogether. Thus, at low energies 
the vector form factors are practically irrelevant. 

We conclude that at GeV energies, the effect of a non- 
zero 77 can be mimicked by a larger value of Ma. On 



the other hand, for E < few 100 MeV, the cross section 
depends only on the G'a(O) and the interpretation of a 
precise measurement is much simpler, leading directly to 
a measurement of 77. These are the reasons of interest 
in considering experiments with low energy neutrinos, 
despite the fact that the neutrino-proton cross section 
diminishes strongly (quadratically) with neutrino energy. 

Observable Spectra in Ultra-pure Scintillators: 
Let us consider an intense beam of neutrinos produced 
from pion decay at rest according to the usual decay 
chains 7r+ — > /i+ + e'^ + h'e + 9^ + v^^. The en- 

ergy spectra of these neutrinos are very well-known and 
include a monochromatic line for the muonic neutrinos at 
(to^ — 77z^)/2TO7r = 29.8 MeV and two continuous spectra 
with energy between < < {rnj^ — m1)/2mf, = 52.8 
MeV. We assume 
the same number of i^e, 



4.9 X 10 pions/yr and thus, 
, and Dp^. Such a pion produc- 
tion rate can be provided, for example, by a cyclotron 
with 800 MeV kinetic energy of protons and with a peak 
power of 1 MW [Tg. 

We calculate the neutral current interactions of these 
Ue, J^/j and with the protons of a ultra-pure scintillator 
detector of mass 1 kton located at a distance of 100 me- 
ters and with the low energy performances as of Borexino 
detector. The expected distribution of the elastic scat- 
tering events, 7^(i^) +p — !■ I'iv) + p, in the kinetic energy 
(4) of the final state proton Tp = E'^ — nip can be written as 



dN 
dTr, 



NpJdE, (^^{E,,Tp) $,(£;,) 



(7) 



+ ^(i?.,Tp) $,(i?.) 



where Np is the number of protons in the scintillator, 
CTjx^p are the cross sections and <I>y and $p are the total 
neutrino and antineutrino fluences (i.e., time integrated 
fiuxes) differential in the neutrino energy E,y. The ki- 
netic energy Tp of the proton in the final state and the 
minimum energy of the neutrino in the initial state are 
related as: 



E„. 



^ (277^p + Tp)Tp 



(8) 



The proton kinetic energy is converted into an 'equiva- 
lent' detectable electron energy Ed, the part of Tp that 
goes into scintillation light. This can be calculated ac- 
cording to the empirical Birks formula, sometimes called 
'quenching' function, [51 IT71ITB] 



Ed 



dT 



1 + k-B (dTp/dx) 



(9) 



where (dTp/dx) is the stopping power of the proton that 
depends on the chemical composition of the detector. 

For the purpose of the analysis we use a Gaussian en- 
ergy resolution a{Ed) = 50 keV-\/£^£i/McV similar to 
Borexino [T2]. We assume the chemical composition 
CqHi2, finding that the number of protons in 1 kton is 
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FIG. 2. Expected event distribution in the detectable energy 
Ed in 1 kton of scintillator as Borexino. The thick solid and 
dotted black lines show the event spectrum for neutral cur- 
rent elastic interactions of neutrinos and antineutrinos with 
protons. The dashed blue line represents the cosmogenic C 
background. The magenta, pink, gray and black solid lines 
represent ^B, CNO, ^ Be and PeP solar neutrino spectrum. 
The purple, red and green dashed lines represent the back- 
grounds due to '^^'^Po, ^^'^Bi and Kr, respectively I21f . 



Np = 6.02 X 10^1 (this increases only by 20% for the 
composition C6H5C„H2„+i with n = 12 considered for 
LENA d^). The Birks Law is obeyed at the percent 
level with a Birks constant of = 0.0175 [20]. There- 
fore, the detection energy threshold E^^ = 0.25 MeV im- 
plies that the equivalent proton kinetic energy threshold, 
calculated with Eq. is about 1.3 MeV. Using Eq. 
we see that we are sensitive only to the highest part of 
the neutrino energy spectra, namely, E^, > 25 MeV. 

The differential spectrum dN/dE^ expected for 1 year 
of exposure is shown in Fig. [2j The signal due to NC vp 
interactions is plotted using three black thick lines. They 
represent the events expected assuming 77 = 0.0,0.12 
(best fit of [3]) and 0.25. From Fig. |2]we see that the 
energy window of interest for the detection of the elastic 
scattering interactions is Ed = (0.25 — 2.2) MeV. The 
number of signal events expected in this energy window 
and for one year of exposure is considerable 
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However, the same energy window will also contain 
several solar as well as radioactive background events. In 
particular solar neutrino events from ^_Be, CNO, PeP 
and are not negligible; we expect about 28.6, 3.0, 
2.1 and 0.5 counts per day (cpd) in each 100 tons of 
scintillator respectively. Moreover, the backgrounds due 
to radioactive sources as ^^'^Po, '^^^Bi and ^^Kr, as well 
as the cosmogenic ^^C have to be considered. Assuming 
the counting rate of the present Borexino configuration 



|21) . we expect that each 100 tons of scintillator yield 
10^ cpd from ^^°Po, 15 cpd from '^^^Bi, and finally, 25 
cpd from both ^^Kr and ^^C. The spectra for these 
backgrounds have also been displayed in Fig. [2j 

Note incidentally that the Borexino detector is still 
improving on the background, and in the most recent 
runs, the ^^'^Po contaminant has already been halved 
j22| . The contribution of the ^^C is particularly criti- 
cal for our purposes, and depends strongly on the depth 
of the underground site of the detector [23]. The value 
of this cosmogenic component is estimated to be 3 times 
smaller in LENA [THj and much less in SNO-I- [H]. How- 
ever, in the present study we will adopt the conservative 
background levels discussed above. Having outlined the 
expectations, we discuss now the sensitivity of the desired 
measurement. 

Sensitivity to proton strangeness r/: From Fig. [2] 
it is evident that background sources can pollute the ex- 
traction of the signal. However, following 25J, we pro- 
ceed to extract the signal from the total number of ob- 
served events, measuring the background events and sub- 
tracting them. 

If we collect N = S -\- B events in a time T and A'o — 
B-Tq/T events in a time Tq, when the signal is absent, we 
expect S — N—Nq-T/Tq signal events. Since both N and 
iVo are subject to Poissonian fluctuations, the fractional 
uncertainty caused by the statistical effects is given by: 
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(10) 



If observations are binned in energy, we define 2/AS'^ = 
d^x^)/dS^ using = {S^ - S,)yASf and S, = Sp,, 
with Pi — 1. In the minimum of the x^. Si — Si, 




-1/2 



(11) 



This is smaller than Eq. ( 10 1 where we set S — J2i ^'^d 
B — J2i-^i' unless the signal and the background have 
the same shape. For constant background-signal ratio, 
Bi/Si =constant, both the results coincide. 

From Fig. [2] it is evident that there are two compara- 
tively clean energy windows to extract the signal, one be- 
tween (0.65 - 1.1) MeV and another between (1.73 - 2.2) 
MeV respectively. Using T = Tq =1 year, we expect 
Si = 28900 and Bi = 37200 events from the first en- 
ergy window and ^2 — 1900 and B2 — 2400 events from 
the latter one. Since we have Bi/Si ~ B2/S2 we can 
use Eq. (10 1 obtaining AS/S = 1.1 %, dominated by 
the data in the first energy window. If we consider that 
for a typical pulsed signal [T3], Tq is 4 times larger than 
T — 1, the error decreases to AS/S = 0.9%; a similar 
improvement is obtained if we measure the background 
over the whole range 0.65 < E^ < 2.5 MeV, increasing 
the statistic and thus reducing its error. 



4 



The NC cross section cr at low energy (and therefore 
the signal event rate S) is dominated by the axial form 
factor, scaling as (1 + 77)^, as is clear from Eq. So we 
can relate the signal sensitivity to the one on 77 through 
SS = Aa/a = 2Ari/{l+ri). This means that the sensitiv- 
ity AS/S = 1.1% on the NC signal implies that proton 
strangeness can be measured with an absolute error of 
A77 = 6 • 10~^ for 7] = 0.12, i.e., a 5% fractional un- 
certainty. For all values of 77 in the allowed experimen- 
tal range, the error is more than 10 times smaller than 
Eq. (|5|. 

Finally, we show how to attain the same sensitivity 
with different experimental parameters. Let us write the 
uncertainty in 77 as follows. 



A.,^i±i' /"'^(;;^^) (12, 

2 ^ e^p S 

Here e — [M' /M){T' /T) rescales the exposure, propor- 
tional to the mass M' and to the time of data taking 
T'; T = {T m)HT/To) rescales the fraction of data tak- 
ing time that includes the signal; ip = {P' / P){D / D'Y 
rescales the flux, accounting for the new power P' and the 
new distance D' . Let us adopt the operational param- 
eters of [13], i.e., assume a distance from the cyclotron 
of D' = 700 m, a power increased by a factor of 5 and 
T' = TqI A = 1 year. The flux (and the signal) decreases 
by one order of magnitude, for it increases with the power 
but decreases as the distance squared. In order to com- 



pensate for this, we need a larger exposure. E.g., with a 
mass of 50 kton as in the LENA proposal [19] we obtain 
again AS/S = 1.1 % after 1 year. Note that the selection 
of the optimal distance between source and detector will 
require to minimize the beam related backgrounds. 

Summary: In this work, we explored the possibil- 
ity to study proton strangeness using low energy neu- 
trinos. We discussed the potential of exploiting a syn- 
ergy between artificial neutrino beams from pion decay 
and ultra-pure scintillating detectors. Incidentally, such 
complex will allow us to quantify precisely the response 
of underground neutrino detectors, to search for sterile 
neutrinos and to investigate many more issues. 

We have shown that, by using reasonable assumptions 
on the experimental parameters, it is possible to identify 
two energy windows that allow to measure many tens of 
thousand of signal events. They imply a statistical error 
on proton strangeness one order of magnitude smaller 
than obtained by BNL 734, Eq. It is important to 
emphasize that the signal we discussed should be col- 
lected above 0.6 MeV, and thus, it does not require the 
extreme low energy threshold reached by Borexino detec- 
tor. 
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